Cosmological Consequences of Slow-Moving Bubbles in First-Order Phase Transitions 
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In cosmological first-order phase transitions, the progress of true- vacuum bubbles is expected to be 
significantly retarded by the interaction between the bubble wall and the hot plasma. We examine 
the evolution and collision of slow-moving true-vacuum bubbles. Our lattice simulations indicate 
that phase oscillations, predicted and observed in systems with a local symmetry and with a global 
symmetry where the bubbles move at speeds less than the speed of light, do not occur inside collisions 
of slow-moving local-symmetry bubbles. We observe almost instantaneous phase equilibration which 
would lead to a decrease in the expected initial defect density, or possibly prevent defects from 
forming at all. We illustrate our findings with an example of defect formation suppressed in slow- 
moving bubbles. Slow-moving bubble walls also prevent the formation of 'extra defects', and in 
the presence of plasma conductivity may lead to an increase in the magnitude of any primordial 
magnetic field formed. 

PACS numbers: 98.80.Cq, 11.27. +d, 64.60.Qb, 64.60.-i 



I. INTRODUCTION 

According to standard cosmology, the early universe 
is expected to have undergone a series of symmetry- 
breaking phase transitions as it expanded and cooled, 
at which topological defects may have formed JlJ . Phase 
transitions are labelled first- or second-order, according 
to whether the position of the vacuum state in field space 
changes discontinuously or continuously as the critical 
temperature is crossed. A first-order phase transition 
proceeds by bubble nucleation and expansion. When at 
least (4— n) of these bubbles collide (for n — 0, 1 or 2), an 
n-dimensional topological defect may form in the region 
between them. 

In recent years there has been considerable interest in 
the formation of defects in first-order phase transitions, 
in particular the validity of the so-called geodesic rule. 
The geodesic rule, first stated by Kibble 0, predicts that 
after a two-bubble collision the phase of the scalar field 
interpolates continuously between the values in each bub- 
ble, along the shortest path in field space. Early anal- 
ysis H confirmed the geodesic rule for defect formation 
in both global and local theories, albeit using a planar 
approximation and neglecting the effect of the surround- 
ing plasma. In later work, the finite conductivity of the 
plasma was considered Ji| for local theories, as was the 
effect of slow-moving (i.e. speeds less than the speed of 
light) bubble walls in theories with a global symmetry 
Q . These analyses confirm defect formation in first-order 
phase transitions, but make conflicting claims about the 
number of defects actually formed. In this paper we in- 
vestigate this issue. Unlike previous work, we include 
the effect of slow-moving bubble walls in both the global 
and local cases. We use our results to make qualitative 
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comparisons between defect densities formed in global 
and local theories, and by slow-moving and fast-moving 
bubble walls. As well as the consequences for defect for- 
mation, we also consider the implications of slow-moving 
walls for the formation of primordial magnetic fields at a 
first-order phase transition. 

We take as our model the simplest spontaneously- 
broken gauge symmetry: the Abelian Higgs model, which 
has a local £7(1) symmetry, with Lagrangian 

L = {D^)\D^) - -F^F» V - V(&$), (1) 



where D^<& 
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The detailed form of the effective potential U(|$|) will 
depend upon the particular particle-physics model being 
considered, but in order to be able to study the generic 
features of a first-order phase transition we shall take V, 
following Ferrera and Melfo j{| , to be 



U($) = A 
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where in a realistic model, e = e(T) cx (T c — T). V 
has a local minimum false-vacuum state at $ = which 
is invariant under the E/(l) symmetry, and global min- 
ima true- vacuum states on the circle |<i>| = pt v (»?/4) (3 + 
e + \/l + 6e + e 2 ) which possess no symmetry. The di- 
mensionless parameter e is responsible for lifting the de- 
generacy between the two sets of minima - the greater 
e, the greater the potential difference between the false- 
and true-vacuum states, and hence the faster the bub- 
bles will accelerate. By making the field and coordinate 
transformations 
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it is possible to set A and 77 to unity, so that the potential 
is parametrized only by e, and hereafter we shall use these 
transformed variables. 

The bubble nucleation rate per unit time per unit vol- 
ume is given by the 'bounce' solution of the Euclidean 
field theory Ignoring quantum fluctuations, the phase 
is constant within each bubble, and uncorrelated be- 
tween spatially-separated bubbles. Any non-zero gauge 
fields in the nucleation configuration will make a contri- 
bution to the action and hence the nucleation of bubbles 
with non-zero gauge fields is exponentially suppressed. 
When three or more bubbles collide, a phase-winding 
of 2im can occur around a point, which by continuity 
must then be at <& = 0. In three spatial dimensions, this 
topologically-stable region of high-energy false vacuum is 
string-like - a cosmic string. 

The formation and evolution of cosmic strings have 
been studied in great detail. Cosmic strings have been 
evoked as, amongst other things, possible seeds for cosmic 
structure formation, sources of cosmic rays, gravitational 
radiation and baryogenesis (see, e.g. (!]]). In order to be 
able to assess the significance of cosmic strings in the 
evolution of the early universe, it is important to be able 
to estimate the initial defect density accurately. This de- 
pends on how the phases between two or more bubbles 
interpolate after collision. In particular, although strings 
are in general formed when three or more bubbles col- 
lide, a simultaneous three-bubble collision is unlikely - 
one would expect in general two-bubble collisions, with 
a third, or fourth bubble colliding some finite time later. 
If the phase inside a two-bubble collision is able to equi- 
librate quickly, and before a third bubble arrives, there 
may be a strong suppression of the initial string density. 
The effect of phase equilibration on the initial defect den- 
sity was first investigated by Melfo and Perivolaropoulos 
JtJ. They found a decrease of less than 10%, in mod- 
els which possess a global symmetry and with bubbles 
moving at the speed of light. 

The above description of defect formation, however, 
ignores any effect that the hot-plasma background may 
have on the evolution of the Higgs field, which may be sig- 
nificant in the early universe. Real-time simulations Iq] 
and analytic calculations |J for the (Standard Model) 
electroweak phase transition predicted that the bubble 
wall would reach a terminal velocity v ter ~ 0.1c. The 
reason for this is simple: outside the bubble, where the 
(SU(2) x U(l)) symmetry remains unbroken, all fields 
coupled to the Higgs are massless, acquiring their mass 
from the vacuum expectation value of the Higgs in the 
spontaneously-broken symmetry phase inside the bub- 
ble. Particles outside the bubble without enough energy 
to become massive inside bounce off of the bubble wall, 
retarding its progress through the plasma. The faster the 
bubble is moving, the greater the momentum transfer in 
each collision, and hence the stronger the retarding force. 
Thus a force proportional to the bubble-wall velocity ap- 
pears in the effective equations of motion. 

Ferrera and Melfo || studied bubble collisions in such 



an environment, for theories which possessed a global 
symmetry, and found that decaying phase oscillations oc- 
cur inside a two-bubble collision, leading to a suppression 
of the defect formation rate [[to| . Kibble and Vilenkin Q| 
studied phase dynamics in collisions of undamped bub- 
bles in models with a local symmetry, and found, analyti- 
cally, a different kind of decaying phase oscillation. When 
the finite conductivity of the plasma was included, these 
oscillations were found not to occur. However, Kibble 
and Vilenkin did not consider the behaviour of the phase 
after collisions of bubbles moving at speeds slower than 
the speed of light. Moreover, because of the symmetry 
assumptions made in their calculations, their results can- 
not be simply extrapolated to the slower-moving case. 

The behaviour of the phase inside bubble collisions in 
local theories where the bubbles move at the speed of 
light, and in global theories with slow-moving bubbles 
has been considered. However, the most realistic sce- 
nario cosmologically - a gauge-theory phase transition 
where the bubbles are slowed significantly by the plasma 
(as might be expected at the electroweak- or GUT-scales) 
- has not been studied. This paper presents the results 
of our investigations into what happens in theories with 
a local symmetry, where the bubbles are moving at ter- 
minal velocities less than the speed of light. Our 3 + 1- 
dimensional simulations indicate that, for slow- moving 
bubbles, phase oscillations of either of the types described 
in H or |U do not occur, before the effect of the plasma 
conductivity is even considered. We therefore expect that 
(a) fewer defects would form in a phase transition where 
the 'Higgs' field is coupled to a gauge field than in a 
global-symmetry phase transition, and (b) in local the- 
ories, fewer defects would form in slow-moving bubbles 
than fast-moving ones. 

We should note in passing that we have ignored the 
effect of the expansion of the universe in our work. This 
is a good approximation for phase transitions which take 
place at late times, like the electroweak phase transition. 
At phase transitions which occur earlier, however, the 
Hubble expansion may have a significant effect on bubble 
and phase dynamics. This topic deserves consideration 
on its own, and work is currently in progress p"l| ] . 

In the following section, we describe the effects of a 
slow-moving bubble wall on phase dynamics inside bub- 
bles collisions in theories which possess a global symme- 
try. In section III, we discuss phase equilibration in theo- 
ries with a local symmetry and present our new results in 
the case of slowly-moving local-symmetry bubbles. Our 
conclusions are supported by examples of defect forma- 
tion suppressed in slow-moving bubbles. We show that 
the 'extra defects' found in |lj] do not occur in heavily- 
damped environments. In section IV we discuss the for- 
mation of a primordial magnetic field. We show that the 
presence of the plasma conductivity results in a larger 
magnetic field for fast-moving bubbles. For slow-moving 
bubbles, the plasma conductivity stops the field dispers- 
ing. A larger magnetic field could also result in this case. 
A discussion of our results and conclusion are presented 
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in section V. 



II. GLOBAL SYMMETRY 

If the gauge coupling e is set to zero, we have a theory 
with a global U(l) symmetry 



(6) 



By writing $ = pe 10 , the equations of motion for the 
modulus p and phase 9 of the Higgs field are 



dV 



p"-(Wp = -— 

d^[p 2 d^] =0. 



dp 



(7) 
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If the potential difference between the true- and false- 
vacuum states is much smaller than the height of the bar- 
rier separating them, the field equations may be solved 
using the 'thin-wall' approximation ||, by setting e = 0. 
For our potential (|2|), this yields 



1*1 = 2 
1 1 2 



1 + tanh 



(a - Rq) 



(9) 



where s 2 = x 2 — t 2 and Rq is the bubble radius on nucle- 
ation. 

Note that 9 = constant trivially satisfies the phase 
equation (j^), and so if the phase is initially constant 
within each bubble, as we shall assume, there are no 
phase dynamics until the bubbles collide. 

As described in the introduction however, we would 
like to investigate the behaviour of the phase in collisions 
of slow- moving bubbles. For a given theory, by consid- 
ering the Boltzmann equations for scattering off of the 
Higgs field, it is possible to calculate the terminal veloc- 
ity of the bubble wall [Q. Since we are not concerned 
here with the parameters of a specific particle-physics 
model, we choose instead to use a single damping pa- 
rameter r to model the interaction of the Higgs with the 
plasma. In the introduction we claimed that the plasma 
would introduce a term proportional to the bubble-wall 
velocity into the equations of motion. Since the phase 9 
of the Higgs field is not affected by the effects described, 
we assume that the plasma couples only to the modulus 
p. We then have effective equations of motion 



p-p" + Tp-{d»9) 2 p 

d»[p 2 d^9] =o. 
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A damping term of this form has been used by several 
authors || , |Q , |Q , and has also been derived from the 
stress-energy of the H igg s, assuming a coupling to the 
plasma |15[. Heckler |14| estimates V ~ gj^Tc for the 
electroweak phase transition, by comparing the energy 



generated by the frictional damping with the pressure on 
the wall due to the damping. 

The effect of this damping term is that instead of ac- 
celerating up to the speed of light, the bubble walls reach 
a terminal velocity v ter < c. By making the ansatz 
p = p [x — xq {t)] , the terminal velocity can be calculated 
Q by integrating the equation of motion for p 
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where Ay is the difference in potential energy between 
the true- and false-vacuum states. Assuming that the 
wall has a Lorentz-contracted, moving profile of the form 

(|) 1 
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1 + tanh 



(r - v ter t - Rq) 



, (13) 



the integral in the denominator of (|T^) can be evaluated 



Expanding 7 = (1 



) 1 ^ 2 , we obtain 
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where A = 6 AV/V\pf v . We have simulated the evo- 
lution of bubbles in such a dissipative environment in 
1 + 1-dimensions. Taking a static profile of the form @ as 
the initial conditions, the terminal velocity of the bubble 
was calculated for a range of values of friction parameter 
r. The accuracy of formula (0), compared with termi- 
nal velocities calculated directly from simulations can be 
seen in Figure 0. This is a very useful result, as from 
it we can dial the input value of T to produce the value 
of v ter corresponding to the particular particle-physics 
model we are interested in. Heckler jl4| , and Ferrera and 
Melfo || obtained a result like ( |l2|) , and Haas jl6| found 
the best-fit equation v ter = A + (1 - A) / (l + BY 162 ) 
from Langevin-equation simulations. However, equation 
(|j"4|), we believe, holds for all of the cases above, provided 
that e is small enough for the 'thin wall' approximation 
to hold, and is more useful when performing simulations. 
For example, it could be applied to the electroweak phase 
transition in the supersymmetric case, if the terminal ve- 
locity of the bubble wall were calculated. 

Ferrera and Melfo j|] described how, in the context 
of a theory with a global symmetry, slow-moving bub- 
ble walls lead to phase oscillations. When two bubbles 
collide, the walls merge. Across the plane (in 3 spatial 
dimensions) of intersection, there exists a phase gradient 
to drive equation (||), and a phase wave propagates into 
each bubble from the centre - see Figures || (b) and || 
(b). As the Goldstone boson is massless, and undamped, 
this wave travels at the speed of light. If the phase dif- 
ference between the bubbles is A9, the phase wave will 
carry a phase difference +A9/2 into one of the bubbles, 
and —A9/2 into the other, equilibrating the phase. If the 
bubble wall is moving at a terminal velocity vt er < c, the 
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FIG. 1. Terminal velocity v te r of bubble walls vs. damping 
coefficient F for e = 0.8 and e = 0.5. The error bars rep- 
resent velocities calculated from simulations, and the circles 
and triangles the values obtained from (|l4|). 

wave will catch up with the bubble wall, and rebound 
- the returning wave will now 'flip' the original phase 
profile. Thus phase oscillations occur inside the merged 
bubbles. Given three or more spatially-separated bub- 
bles whose distribution of phases one would expect to 
generate a vortex on collision, a vortex, an anti-vortex, 
or none at all may form, depending on the profile of the 
phase inside the two bubbles at the moment of collision 
of the third - an example of how phase dynamics can 
affect the defect-formation process. The oscillations are 
damped, because the bubble walls continue to expand, 
increasing the volume over which the finite-energy wave 
must sweep, thus diluting the phase difference carried 
by the wave. Thus the converse of the above statement 
is not true - an initial distribution of phases which one 
would not expect to form a defect, will not produce one 
as a result of phase oscillations. Statistical simulations 
in two dimensions [[lOf have shown that this leads to a 
suppression in the defect-formation rate - the slower the 
bubble walls, the fewer defects are formed per nucleated 
bubble. 



III. LOCAL SYMMETRY 

A. Phase dynamics inside two-bubble collisions 

Including the gauge fields in our model, the field equa- 
tions become 



p-p" - {3^9 - eA^fp 



dV 
"dp 



0"[p 2 (0 M 0- eJ 4 M )] =0 
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Since we now have a local U(l) symmetry, the phase 
9 can be arbitrarily re-defined at any point in time by a 
gauge transformation, and so we need a gauge-invariant 
notion of phase. We define, following Kibble and Vilcnkin 
the gauge- invariant phase difference between two 
points A and B 



A6 = 



dx l (di — ieAi) . 



(18) 



where i — 1,2,3 and the integral is taken, for simplicity, 
along the straight line joining A and B. 

For bubbles which move at approximately the speed of 
light, it is possible to greatly simplify the field equations. 
If we consider a two-bubble collision, in a frame where the 
bubbles are nucleated simultaneously, by assuming that 
the bubbles instantly propagate at the speed of light, it 
is possible to impose 5*0(1,2) Lorentz symmetry on the 
field equations. Thus the fields are functions of z and 
t 2 = t 2 — x 2 — y 2 only. With this assumption, and a step- 
function ansatz for the phase 9 at the time of collision, 
Kibble and Vilenkin j| solved the field equations for A9 

2R ( 1 \ 

AO = — O coserjCt- R) H sin er] (t - R) . (19) 

t \ er\R J 

where 29$ is the initial phase difference between the 
spatially-separated bubbles and R is their radius on col- 
lision at t = 0. 

Equation ( [l9| ) describes decaying phase oscillations, 
the time scale of equilibration determined by the initial 
phase difference and radius of the colliding bubbles, the 
frequency of oscillation by the gauge-boson mass. These 
oscillations, and the accuracy of this formula for small 
initial phase differences, have recently been confirmed in 
simulations fl7|| . 

However, the assumption that the bubbles move at, or 
close to the speed of light, does not appear to be realistic 
Q. In this case, the symmetry assumptions made in Q] 
are no longer valid. Moreover, it is not possible to replace 
the coordinate r 2 = t 2 

T' 2 = (Vtert) 2 - £ 

the Higgs field, the bubble wall, is constrained in this way 
- the phase and gauge fields are still free to propagate 
causally. 

In order to investigate whether phase oscillations - 
which occur in the global theory with slow-moving bub- 
bles, and in the local theory with fast-moving bubbles - 
still occur in the local theory when the bubbles expand 
slower than the speed of light, we include the dissipa- 
tion term Tp into the equation for the modulus of the 
Higgs field p, without coupling it to the phase or the 
gauge fields. This is motivated in the same manner as 
described in the global case. A term proportional to p is, 
of course, (7(1) gauge-invariant. 

Since there is no longer an obvious simplification of the 
equations of motion which might lead to an analytic so- 
lution, we turn to computer simulations. The equations 
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FIG. 2. Undamped bubbles passing through one another: a) with no phase difference, and b) with phase difference A# = 7r/3. 
In both cases £ = 0.1. The direction of the arrows gives the phase # and their length is proportional to the modulus p. Note in 
case b) that a phase wave propagates at the speed of light from the point of collision into the bubble interiors. In this case, the 
symmetry-restored region is smaller than in a), because the phase wave is able to carry away energy from the collision point, 
which is not possible in case a), where there is no phase difference between the two bubbles. 



of motion were discretized in the gauge-invariant way de- 
scribed in Jl8fl , choosing the temporal gauge Ao = in 
order to make the time evolution trivial. We used a lat- 
tice of size 200 3 and a lattice spacing a = 0.5 - tests were 
performed on lattices with spacing down to a = 0.1 giv- 
ing no qualitatively-different results. The time evolution 
was performed using a fourth-order Runge-Kutta algo- 
rithm. We took as initial conditions a static profile of 
the form (^|) for p = |$|, for two bubbles of radius R = 5, 
with phases 9 = and 9 = 2ir/3, centred at (±8,0,0). 
We choose to ignore any primordial magnetic field and, 
since the nucleation process is not expected to generate 
non-zero gauge fields (see Introduction), set all the gauge 
fields to zero initially. 

The results of the simulations are displayed in Figures 
|], ^| and For the sake of clarity and to aid comparison 
between the different cases, we have chosen to present our 
results in terms of the evolution with time of the gauge- 
invariant phase difference A9. We evaluated A9 between 
the centres of the two bubbles, though the qualitative 
behaviour was found not to change when it was calculated 
between different points. 

Figure |J (a) shows the behaviour of the gauge- invariant 
phase difference for bubbles moving at the speed of light 
- the decaying oscillations calculated by Kibble and 
Vilenkin in the local case. In the global case, e = 0, 
we find that the phase does equilibrate, but on a much 
longer time-scale. Thus we would expect that for fast- 
moving bubbles, fewer defects are formed in local theories 
than global ones, since in order to form a defect a phase 
difference inside the two merged bubbles must be present 
when a third bubble collides. 

In Figure ^ (b) we plot A9 for slower-moving bub- 
bles. For e = 0, we confirm in 3 + 1-dimensions the de- 



caying phase oscillations described by Ferrera and Melfo 
Q and observed by them in 2 + 1-dimensions. These 
oscillations are killed by adding in gauge fields - for a 
fixed bubble-wall velocity, the stronger the gauge cou- 
pling, the less time the gauge-invariant phase difference 
is non-zero, and hence the less likely a third collision will 
occur in time for a defect to form. Thus we would ex- 
pect a lower defect-formation rate in local theories with 
slower-moving bubble walls. 

Figure |E| illustrates our findings - it shows a cross- 
section through a non-simultaneous three-bubble col- 
lision, after all three bubbles have merged. In each 
case, the bubbles of initial radius R = 5, centred at 
(±8, 0, -10) and (0, 0, 10), were given phases 9 = -tt/2, 
and 27r/3. For identical initial conditions, we see that in 
the fast-moving case a vortex is formed, but when the 
bubbles are slowed down, the phase difference between 
the two bubbles has equilibrated by the time the third 
bubble collides, and no defect is formed. 

In any cosmological phase transition where the bubble 
wall is significantly slowed down, we may also expect the 
plasma to have non-zero conductivity, which will affect 
the evolution of the fields and so needs to be considered in 
any attempt at a realistic model. We have simulated the 
effects of the finite-conductivity of the plasma, by adding 
on to the right-hand side of the gauge field equations ( |l7| ) 
a conduction current jg, whose spatial part is given by 

jc = crE. (20) 

The corresponding charge density p c is fixed by the con- 
tinuity relation d^j£ = 0. For large values of the con- 
ductivity, it has been shown jl| that the oscillations in 
the gauge-invariant phase difference, which took place 
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FIG. 3. Damped bubbles bubbles merging: a) with no phase difference, and b) with phase difference AS = tt/3. In both 
cases £ = 0.1. Note in case b) the phase wave propagating inside the bubbles at the speed of light after collision, faster than 
the bubble walls. This phase wave will catch up with the outer wall of the bubble (not shown here), and rebound, producing 
phase oscillations. 



in fast-moving bubbles with cr = 0, are exponentially 
damped. 

Figure ^| (a) shows the evolution of the gauge- invariant 
phase difference in this case, where the walls are moving 
at the speed of light, for three different values of the con- 
ductivity cr. We confirm that, as a increases, the phase 
oscillations are more heavily suppressed, with practically 
no oscillations occurring for a > 0.5. In Figure || (b), 
we present the results of our simulations for slow-moving 
bubbles. For a — 0, we have the case considered above 
- heavily suppressed oscillations. Increasing a merely 
serves to increase the suppression of phase oscillations: 
no new effect is observed. 

Whereas it is already known that for bubbles moving 
at the speed of light, phase oscillations can be killed by 
a high conductivity || , it is clear from our work that in 
slower-moving bubbles, the same effect can be obtained 
by a much lower value of a. 

B. Extra defect production 

An interesting consequence of slower-moving bubbles 
concerns the issue of 'extra' defect production at a two- 
bubble collision. Hawking, Moss and Stewart (Til first de- 
scribed (by energy considerations) how two true-vacuum 
bubbles travelling at nearly the speed of light would 'pass 
through each other', leading to the temporary restora- 
tion of the spontaneously-broken symmetry in a region 
between the two bubbles. This is illustrated in Figure 
U - the two bubbles collide, and bounce off of (or pass 
through) each other, producing a region of <& = false 
vacuum inside the merged bubbles, which decays via os- 
cillations of the bubble walls into the true- vacuum state. 
The size of the symmetry-restored region, and the time 



taken to decay completely to the true vacuum depends on 
the initial phase difference between the two bubbles, and 
the asymmetry parameter e. Copeland and Saffin (T^ ] 
showed how this could lead to the formation of 'extra' 
- in the sense that a defect would not be expected from 
the initial distribution of phases - flux-tube vortices in 
a gauge theory, around these regions of temporarily re- 
stored symmetry, and hence to an increase in the initial 
defect density after a phase transition. 

Our simulations show - see Figure || - that dissipa- 
tion prevents this bouncing, or passing-through, of the 
bubbles. The excess energy, which would cause the sym- 
metry restoration, is dissipated away by the plasma, and 
the bubbles simply merge. Thus there is no symmetry- 
restored region around which a non-zero winding of the 
phase can occur, and so no 'extra' defects would be 
formed. 



IV. MAGNETIC FIELDS 

Another consequence of first-order phase transitions 
which may be cosmologically significant is the generation 
of primordial magnetic fields. Galaxies are observed to 
have magnetic fields B ga \ ~ 10~ 6 G, coherent over large 
scales. Given a small initial seed field a dynamo mecha- 
nism, powered by the differential rotation of the galaxy in 
combination with the small-scale turbulent motion of the 
ionized gas, could generate the observed galactic fields. 
Many mechanisms for producing such a seed field have 
been proposed, one being bubble collisions at a first-order 
phase transition (see e.g. [^0| and references therein). It 
had been believed that it was not possible to generate a 
seed field of sufficient magnitude at the electroweak phase 
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FIG. 4. Gauge-invariant phase difference between a) two bubbfes moving at the speed of light , F = 0, and b) two slow-moving 
bubbles, T = 2. 



transition for the dynamo mechanism to explain galac- 
tic magnetic fields as large as 10 _6 G. However a recent 
paper by Davis, Lilley and Tornkvist [pl| showed how, 
in a universe with a low matter density and in particu- 
lar a positive cosmological constant, a dynamo mecha- 
nism may be able to generate observed galactic magnetic 
fields from a much smaller seed field. As a consequence, 
electroweak-scale magnetic fields may be viable primor- 
dial seed fields, and it is of interest to consider the effect 
of slow-moving bubble walls and finite plasma conductiv- 
ity on the generation of magnetic fields. 

If the gauge fields are set to zero initially, it can be 
seen from the equations of motion (1^), ( |l6| ) and (17), 
that non-zero gauge fields can only be generated where 
there exist spatial phase gradients, that is after the col- 
lision of two or more bubbles. After the collision of two 
bubbles, a loop of magnetic flux is generated around the 
circle of intersection [Q . The amount of flux generated 
is given by the integral of the gauge field A around any 
loop which passes outside the bubbles. This is the same 
for all bubbles, regardless of size or speed 



26n 



(21) 



and our simulations confirm this. When a third bubble 
collides, the fluxes combine, and if there is a phase wind- 
ing of 27T around the centre, one flux quantum 2ir/e will 
be trapped. 

If there is no plasma conductivity, the magnetic flux 
generated is free to propagate at the speed of light away 
from the bubble collision. If the bubbles are expanding at 
the speed of light, then the fields can disperse no further 
outwards than the intersection of the bubbles. Copeland, 
Saffin and Tornkvist j p2fl , [jlTj demonstrated how in this 
case two tubes of flux are produced - a 'primary flux tube' 
at the intersection of the collided bubbles and a smaller, 
'secondary' peak of opposite direction following behind 



it - see Figure ^ (a). Our simulations show that when 
non-zero conductivity is included, this secondary peak 
does not occur, and all the magnetic flux is concentrated 
into the primary peak, which is consequently larger - 
Figure (b). We might expect that in this case, a larger 
magnetic field would form as all the flux generated by 
the two-bubble collision is aligned. 

If, however, the bubbles are moving at speeds less than 
the speed of light, the flux is able to disperse into the 
plasma. Unless the bubble nucleation rate is extremely 
high (when a third bubble might be expected to collide 
quickly after the initial collision), no magnetic field will 
be able to form. For slow-moving bubbles, large-enough 
conductivity prevents the flux from dispersing, freezing 
it in to the plasma. Figure ^ (c) shows the magnetic 
field strength formed after a collision of two slow-moving 
bubbles, for a — 5 (in reference an estimate of a = 
T/e 2 is given. The temperature at a phase transition is 
typically T ~ 77, and so a = 5 may be realistic) . It can be 
seen that for this value of a, a magnetic field does form, 
but it is spread through the inside of the bubble, rather 
than being concentrated in one or two narrow peaks at 
the wall, as seen in the fast-moving case. The height 
of the peak is lower by an order of magnitude - this is 
finite, rather than infinite conductivity and so some flux 
still escapes. 

We note in passing here that provided that the plasma 
dynamics which slow down the bubble walls do not af- 
fect the bubble nucleation rate, the average number of 
bubbles nucleated per unit volume by the time the phase 
transition is completed will increase as the bubble wall 
velocity decreases. That is the average bubble radius on 
collision, the correlation length of the Higgs field £, de- 
creases as the wall velocity decreases. Since the amount 
of flux generated at each collision is independent of the 
bubble radius, slow-moving bubbles will generate more 
flux, and hence a larger magnetic field when coarse- 
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FIG. 5. Phase plot and bubble walls after three-bubble collisions, with phases (bottom left), 27r/3 (top left) and —n/2 
(right): a) with V = a vortex is formed at the centre, and b) with identical initial conditions, but V = 0.5 there is no vortex. 



grained over many bubble radii. 

V. CONCLUSION 

In this paper, we have examined the behaviour of col- 
liding true-vacuum bubbles at a first-order cosmological 
phase transition. In the Abelian Higgs model, strings 
may form at the collision of three or more bubbles, but 
since a simultaneous three-bubble collision is very un- 
likely, the dynamics of the phase inside two-bubble colli- 
sions is crucial - if phase differences between two bubbles 
can be equilibrated quickly, and before the arrival of a 
third bubble, a topological defect will not form. 

The most relevant phase transitions to cosmology in- 
volve gauge fields coupled to the symmetry-breaking 
field. In such phase transitions, the speed of the bub- 
ble walls will be considerably less than the speed of light, 
and yet the phase dynamics of slow-moving bubbles in 
a gauge field have not been considered previously. We 
have thus paid particular attention to the evolution of 
the phase inside collisions of bubbles moving at speeds 
much lower than the speed of light, in a £7(1) gauge the- 
ory. 

In the simplest model, with no gauge fields and where 
the bubble walls accelerate up to the speed of light, the 
phase difference between two points is found to equili- 
brate. In models with a global symmetry where the bub- 
ble walls move slowly, and models with a local symmetry 
where the walls move at the speed of light, decaying phase 
oscillations have been observed. We find that in a U(l) 
gauge theory, with slow-moving bubble walls, these oscil- 
lations are suppressed. On collision of two bubbles, in- 
stantaneous phase equilibration is observed. This would 
lead to a decrease in the initial expected defect density 
compared to the other cases. We have illustrated our 
claims by demonstrating an example of the suppression 



of defect formation in a local theory, due to nontrivial 
phase dynamics. 

When two bubbles collide and merge, there will exist 
phase gradients across the intersection, a potential differ- 
ence. In local theories it is necessary to define a gauge- 
invariant notion of the phase difference, which involves 
the gauge fields. Thus in local theories, the phase dif- 
ference may equilibrate through the generation of gauge 
fields - there is in effect an 'extra channel' for the de- 
cay of the potential difference created on collision. This 
explains why we would expect fewer defects in local the- 
ories than in global ones. In a local theory, the phase 
difference between two bubbles is observed to equilibrate 
more quickly in slower-moving bubbles than in bubbles 
moving at the speed of light. This is due to the fact in 
slow-moving bubbles the rate of generation of phase gra- 
dients is lower, yet the gauge fields are not restricted to 
propagate at the speed of the bubble wall and are thus 
able to equilibrate the phase difference more rapidly. 

If phase equilibration and hence the suppression of de- 
fect formation is aided by the coupling of gauge fields to 
the Higgs, it is interesting to ask whether another scalar 
field x could have the same effect. In order to be able to 
dissipate the potential energy in the phase gradient such 
a field would need to couple to the phase of the Higgs 
field, but also preserve the U(l) symmetry of the La- 
grangian. This can only be achieved (with terms at most 
quadratic in <E> and its derivatives) by Higgs couplings 
proportional to d„ X [^d^ - (5^)$] = d, lX [2ip 2 d li ei 
but in this case x is effectively a gauge field 0. It is 
possible that fermion couplings to the Higgs field would 
aid phase equilibration, but unfortunately this cannot be 
simulated easily. 

We predict that fewer defects will form in gauge the- 
ories than global-symmetry theories since the phase dif- 
ference is non-zero for less time after collision in gauge 
theories. In fact, if the bubble nucleation rate is low 
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FIG. 6. Gauge-invariant phase difference between two bubbles moving at the speed of light, with conductivity a = 0, 0.1 and 
0.5. In each case, the initial phase difference is and e = 0.5. In Figure a) T — 0, and in Figure b) F = 0.5. 



enough, it might be possible to effectively rule out the 
formation of defects, solely on phase-dynamical grounds, 
though percolation could presumably still be achieved. 
This is a very interesting prospect, which could have sig- 
nificant implications for cosmology - it may be possible 
for example to circumvent the monopole problem with- 
out needing inflation if defect formation is dynamically 
suppressed in this way. 

It has also been seen how it is unlikely that 'extra 
defects', caused by bubbles bouncing off of each other 
on collision, will be formed in cosmological phase transi- 
tions, since the bubbles are retarded sufficiently by the 
plasma for no such bouncing to occur. First-order phase 
transitions can also generate a primordial magnetic field, 
which may seed the galactic dynamo and hence be re- 
sponsible for the galactic magnetic fields observed today. 
A simple qualitative analysis suggests that in fast-moving 
bubble walls, high conductivity (as would be expected 
in the early universe) would lead to the generation of 
a larger magnetic field. Where the bubble walls move 
slower, we have demonstrated that a magnetic field can 
form if the plasma has non-zero conductivity. In this 
case, the smaller average bubble radius on collision may 
cause more flux to be generated, producing a larger mag- 
netic field. This may be significant in helping to beat the 
lower-bound required by the dynamo model in order to 
produce observed fields. 

We have shown qualitatively how we expect the defect- 
formation probability to be decreased by phase equilibra- 
tion in two-bubble collisions. It would be interesting to 
perform a statistical simulation of the type done in , 
but for gauge-theory phase transitions, to see quantita- 
tively how the defect density is affected by the terminal 
velocity of the walls or the introduction of gauge fields. 
We note that the argument given in the Introduction for 
nucleating bubbles with zero gauge fields does not apply 
if there already exists a primordial magnetic field before 



bubble nucleation. In this case, it is not at all clear what 
the preferred nucleation field configuration would be, and 
we believe that a study of bubble nucleation in the pres- 
ence of a magnetic field would be worth while. It is also 
of some interest to consider the effect on phase dynamics 
and defect formation of the Hubble expansion, since this 
acts as a dissipation term on the phase as well as on the 
bubble walls. We conclude, though, with a summary of 
our findings. 

In gauge theories, more defects are formed by fast- 
moving bubble walls than by slower ones. In global the- 
ories, the same is true. 

For fast-moving bubble walls, more defects are formed 
in global theories than in local ones. For slow-moving 
bubble walls, the same is true. 
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